Zero-energy states in rotating trapped Bose-Einstein condensates 
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We solve the quasiparticle excitation spectrum of a rotating three-dimensional Bose-Einstein 
condensate whose ground state is pierced by seven singly quantized vortex filaments. We find the 
elementary excitation spectrum of this system to harbour a quasiparticle with zero excitation energy. 
This quasiparticle is degenerate with the rotating ground state and may be described as a massless 
Nambu-Goldstone mode or a Majorana boson. Macroscopic population in such a state corresponds 
to a Bose-Einstein condensate of Tkachenko vortex wave quanta. 
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The significance of zero-energy excitations in quantum 
field theories of both elementary and quasiparticles has 
substantially increased since Majorana predicted the ap- 
pearance of a mysterious massless particle, which acts as 
its own anti-particle pQ. Neutrinos have been suggested 
as candidates for Majorana fermions and more recently, 
non- elementary Majorana quasiparticles have sparked 
major interest due to their non-Abelian exchange statis- 
tics, which could be harnessed for fault tolerant topolog- 
ical quantum computing [2r|4]. Majorana zero modes are 
also predicted to become trapped inside vortex cores in 
chiral p-wave paired superconductors and protected by 
coupling to edge modes [3HS]- Signatures of Majorana 
modes in superconductor-semiconductor nanowires have 
recently been observed [7]. 

Zero-energy excitations have been central to discus- 
sions of field theories with spontaneously broken contin- 
uous symmetries [SHU] ■ These massless modes are called 
Nambu-Goldstone bosons, the photon being a famous 
example. In non-relativistic condensed matter systems, 
phase transitions involving spontaneous symmetry break- 
ing, such as the Bose-Einstein condensation phase tran- 
sition emerge with associated Nambu-Goldstone bosons. 
In superconductors such zero mode becomes gapped due 
to the Coulomb interaction [T3HT5] . Furthermore, the 
mass of elementary particles is thought to emerge due 
to shifts in the energies of zero modes [15--2T3]. Also in 
systems where continuous symmetries are broken explic- 
itly rather than spontaneously, the Nambu-Goldstone 
bosons typically acquire small mass correction and are 
called pseudo-Nambu-Goldstone bosons [21]. 

Here we study a quantum degenerate Bose gas, which 
is exposed to an external vector potential that causes the 
atoms to rotate at an angular velocity f2 around the cho- 
sen z-axis of the three-dimensional harmonic potential. 
By applying the microscopic Bogoliubov-deGennes field 
theory, we show that this symmetry broken system is ca- 
pable of hosting a zero mode whose physical appearance 
is that of a Tkachenko vortex wave. Such an excitation 
mode is degenerate with the rotating ground state and 
emerges as a consequence of the broken continuous rota- 
tional symmetry in the vortex lattice ground state. Since 



there is no energy cost associated with populating this 
excitation mode, we expect these quasiparticles to be- 
come macroscopically occupied forming a Bose-Einstein 
condensate of Tkachenko wave quanta. 

By using the method of second quantization, Bogoli- 
ubov provided an explanation of the phenomenon of su- 
perfluidity in terms of a gas of quasiparticles |22j . In 
this picture, the quasiparticle creation operator ijj q and 
annihilation operator ^/>J obey the equations of motion 

ihdtljjq = Eq^q 
iHdtj>\ = -Eqi>l (1) 

Quasiparticles with zero energy E q = satisfy the rela- 
tion 

fa = i>l (2) 

which is the defining property of a Majorana operator — 
a (quasi) particle creation operator is equivalent to its 
annihilation operator. The peculiar character of Majo- 
rana quasiparticles is not restricted to systems composed 
of fermionic particles. The Majorana zero mode itself is 
neither a boson or a fermion as it satisfies neither com- 
mutation or anticommutation relations and is therefore 
described as an anyon |23j . 

The quasiparticle energies E q can be explicitly cal- 
culated for a given physical system by solving the 
Bogoliubov-deGennes equations [HI [25] 

£(r)tt g (r) - M{r)v q {r) = (E q + n)u q {r) 
C*(r)v q (r) - e i8 M*(r)u q (r) = -(E q - (i)v q (r) (3) 

where £(r) and A4(r) are model dependent operators, 
u q (r) and v q (r) are the quasiparticle amplitudes, and \x 
is the chemical potential. For the boson system stud- 
ied here = 0, whereas for fermion systems the anti- 
commutation relations yield = tt. For systems whose 
finite-energy excitations are fermionic the excitation en- 
ergies E q are measured with respect to the Fermi energy 
E-p, whereas for boson systems they are measured with 
respect to the chemical potential fj, of the condensate. 

Considering a Bose-Einstein condensed system, a 
quadratic approximation to the full many-body Hamilto- 
nian, expressed in terms of the atom field operator ^(r), 
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FIG. 1. Nambu-Goldstone boson cf>(r) of a condensed Bose 
gas rotating at an angular frequency Q = 0.5Suj±. (a) Isosur- 
face of the condensate density |0(r)| 2 with seven quantized 
vortices piercing the condensate forming a sixfold symmet- 
ric triangular array, viewed along the rotation axis, (b) The 
corresponding phase map S(r) = arg(0(r)) shows the seven 
phase singularities at the vortex cores, marked with filled cir- 
cles. 

can be diagonalized by the canonical Bogoliubov trans- 
formation to a quasiparticle basis [22] 

*(r)-0(r) = ^ Ug (r)^- t ;*(r)^ (4) 

where <p(v) is the macroscopic wavefunction describ- 
ing the atoms in a Bose-Einstein condensate. The 
quasiparticle creation and annihilation operators tp q and 
$J preserve the canonical boson commutation relations 

[i'qti'p] = $qp> [i> q Ap] = Wq^p[ = and tnc nnite en- 
ergy quasiparticle amplitudes satisfy the orthonormaliza- 
tion condition [U] J(u*(r)u p (r) - v*(r)v p (r))d 3 r = S gp 
and J (u q (r)vp(r) - u p (r)v q (r))d 3 r = 0. 

We consider a Bose-Einstein condensate of N = 3 x 10 
atoms at zero temperature limit. The 87 Rb atoms of mass 
m are confined in a three-dimensional harmonic oscillator 
potential V(r) = m^^r 2 ^ + lj 2 z 2 )/2, which rotates at an 
angular frequency Q = 0.58wj_ unless otherwise stated. 
This rotation frequency stabilizes the seven vortex array 
as the ground state of this system. The interaction cou- 
pling constant g — 4irh 2 a/m, where a is the s-wave scat- 
tering length of the atoms and the oscillator length ao = 
We have chosen wj_ = lo z /^/S = 2ir x 100 Hz. 
For this system, the operators in Eq. ^ are M(r) = 
# 2 (r) and £(r) = - ^+V(r)+2 5 |<£(r)| 2 -QL z , where 
L z denotes the projection of the orbital angular momen- 
tum operator onto the axis of rotation. The condensate 
order parameter </>(r) satisfies the Gross-Pitaevskii equa- 
tion [261 [27] 

(£(r)- ff |0(r)| 2 - A1 )^(r)=O. (5) 

To draw an analogy with charged superconductors, the 
operator £(r) may alternatively be expressed in terms 
of a vector potential An = m 2 Q,(—ye x + xe y )/h as 
£(r) = [(ifiV - nA n ) 2 + m 2 (u\ - Sl 2 )r\ + m 2 u 2 z 2 + 
4mg\4>(r)\ 2 ]/2m, where k — h/m and h is Planck's con- 
stant. 



FIG. 2. Majorana boson i/>(r) of a condensed Bose gas rotat- 
ing at an angular frequency Q, — 0.58o;x. (a) Isosurface of 
the density |i/>(r)| 2 viewed along the rotation axis, (b) Rela- 
tive phase map M(r) = arg(^(r)) — arg(^(r)) shows 18 phase 
singularities. The six outermost of these are ghost vortices of 
ip(r) and reside at the edge of the condensate. The remaining 
12 phase defects correspond to 6 vortices in the condensate 
ground state <f>(r) and 6 antivortices in the state ^(r). Both 
4>{t) and ip{r) have a vortex in the center of the trap. These 
vortices are of opposite circulation and cancel each other in 
the relative phase map. 

The coupled equations ^ can be shown to possess 
a particle-hole symmetry such that given a solution 
u q (r), v q (r), E q a second solution vl (r), v' (r), E' always 
exists and is obtained by a simple transformation u' q (r) = 
V*(r), v' q (r) = u*(r) and E' — —E*. For fermion systems 
both hole and particle type excitations are treated on an 
equal footing whereas for bosonic systems half of all ex- 
citations have a negative norm and in order to avoid de- 
scribing them as fermionic excitations, they are typically 
discarded. In addition, for a Bose-Einstein condensed 
system, there exists a trivial solution, E\ = to Eq. (J3j|, 
which is obtained by setting U\{t) — v\(r) = </>(r). 
This is the Nambu-Goldstone boson corresponding to 
the spontaneous U(l) gauge symmetry breaking due to 
the formation of the Bose-Einstein condensate. Such a 
zero mode does not obey the bosonic quasiparticle nor- 
malization condition and instead it satisfies the relation 
/ ul(r)ui(r)d 3 r — J Vi(r)vi(r)d 3 r. In addition to the 
trivial solution, other zero modes degenerate with the 
ground state and which satisfy £(r)u q (r) = A4(r)u*(r), 
may exist in the spectrum. 

We have calculated the quasiparticle excitation spec- 
tra by numerically solving Eq. ^ using a parallelized 
Arnoldi iteration method. The three dimensional oper- 
ators are discretized using a finite-element discrete vari- 
able representation [28] to yield a sparse matrix repre- 
sentation for the Bogoliubov-de Gennes operator. Before 
solving the Bogoliubov-de Gennes problem, the conden- 
sate ground state <p(r) and the corresponding chemical 
potential fi are obtained by solving the Gross-Pitaevskii 
equation ([5| using an over-relaxation method. The self- 
consistency of the calculated spectrum is guaranteed by 
the fact that we obtain the E\ = Nambu-Goldstone 
solution with a numerical precision \Ei\ < 10~ 9 huij_. 

Figure [Ha) shows a condensate density |</>(r)| 2 iso- 
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surface and (b) the corresponding phase map S(r) 
for the calculated ground state wavefunction </>(r) = 
\4>(r)\ exp(iS(r)). The chemical potential of this conden- 
sate is [i = (<f>\£(r) — g\(p(r)\ 2 \<f)) = 12,0fkj± and the 
orbital angular momentum L = (<f>\ — iH(xd y — yd x )\4>) = 
4.2 hN. The condensate is pierced by 7 singly quantized 
vortex filaments. The circulation § v s (r) ■ dl — k£ of the 
superfluid velocity v s (r) — fi/mVS(r) around the vortex 
cores is quantized in integer t multiples of n = h/m. In 
the ground state of the rotating potential, these vortices 
are arranged in a triangular lattice structure with sixfold 
discrete rotational symmetry. The phase singularities lo- 
cated inside the vortex cores are marked with filled circles 
in Fig. Qb). 

Figure [3] shows the Bogoliubov excitation energies cal- 
culated for the rotating condensate ground state shown 
in Fig. [TJ Energies of the lowest collective quasipar- 
ticle modes are shown as a function of an integer in- 
dex q, which labels the quasiparticles in increasing or- 
der of their energy. The most prominent feature of 
this excitation spectrum is the presence of two zero en- 
ergy modes \Ei^/fu^i_\ < 1CT 9 labelled by q = 1 and 
q = 2. We note that all shown states are duplicated 
in the numerically calculated spectrum due to the afore- 
mentioned particle-hole symmetry of the Bogoliubov-de 
Gennes equations and thus we find 4 zero modes, only 
half of which represent distinct quasiparticles and which 
are shown in Fig. [3] The first one of these modes q = 1 
is the Nambu-Goldstone boson <^(r) always present in 
these Bose-Einstein condensed systems irrespective of 
the value of the vector potential Aq. It corresponds to 
the Bose-Einstein condensate of atoms and is equivalent 
to the state shown in Fig. [I] The presence of this mode 
reflects the spontaneously broken U(l) gauge symmetry 
of the ground state as the phase coherent Bose-Einstein 
condensate is formed. 

The observation of the second zero energy mode q = 2 
corresponds to a massless Nambu-Goldstone mode ^>(r) 
or a Majorana boson. In the condensate ground state 
the continuous rotation symmetry SO (2) has also been 
broken and only the reduced discrete sixfold rotational 
symmetry of the vortex lattice remains. The source of 
this symmetry breaking can be traced back to the vector 
potential Aq in the Hamiltonian. Therefore one might a 
priori expect to find a zero energy boson mode in the ex- 
citation spectrum due to Goldstone's theorem. Instead, 
the zero mode has been observed to acquire a small mass 
correction [29]. Furthermore, this zero mode is different 
from the high angular momentum surface mode, which 
approaches zero energy when new vortices are nucleated 
in the system J3UJ 125 ■ 

The physical character of the quasiparticle state "0(r) 
is that of a Tkachenko vortex wave 32 . It comprises 
of elliptically polarized collective motion of the vor- 
tices around their equilibrium locations [3"2T(4"0"] . This 
mode has transverse and longitudinal components cor- 
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FIG. 3. Quasiparticle excitation energies as a function of an 
integer index q which labels the modes in the order of increas- 
ing energy. The two zero modes are singled out in the insets. 
The energy of both of these excitations is \Ei l 2/foui±\ < 10~ 9 . 
The center-of-mass Kohn mode is marked with a triangle. 
The 7 modes whose energies lie between the condensate q = 1 
and the Kohn mode q — 8 are the vortex displacement modes 
corresponding to the 7 vortex degrees of freedom. 
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FIG. 4. Quasiparticle probability density of the zero energy 
Majorana boson mode ip(r) in the z — plane. The six peaks 
are localized in the vortex cores of the ground state 4>(r). At 
the origin, there is a phase singularity in both <p(r) and ip(r) 
where the probability density vanishes. 

responding to azimuthal and radial motion of the vor- 
tices. The Tkachenko modes of vortex lattices in har- 
monically trapped Bose-Einstein condensates have been 
calculated analytically within a vortex continuum ap- 
proximation and the mode frequencies have been found 
to approach zero in the centrifugal limit Q — > loj_ as the 
system approaches the quantum-Hall regime [40 . This 
is in agreement with the experimental observations for 
large vortex lattices [2PJ. However, the Tkachenko mode 
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has not been predicted to appear at zero energy for sys- 
tems with finite £1 < Moreover, previous numerical 
calculations for purely two-dimensional systems [4TH44] 
or three-dimensional prolate systems [45] have not pre- 
dicted a zero mode other than the condensate Nambu- 
Goldstone boson in these systems. We have performed 
calculations for in the range 0-0.8 for the present 

three-dimensional oblate system and have found a zero 
energy vortex wave only for the case of seven vortices. 
In this context, it is amusing that seven classical vortices 
in a ring arrangement is a marginally stable configura- 
tion having a normal mode with precisely zero frequency 
[461147]. 

We may understand the finite energy of a typical 
Tkachenko mode t^(v) in terms of a coupling to the 
condensate mean-field 0(r). If ip{v) absorbs a Nambu- 
Goldstone boson from (j>(r) it generically changes the or- 
bital angular momentum of the condensate, which leads 
to a changing vortex density and is tantamount to ra- 
dial vortex motion. The restoring force provided by the 
confinement results in the sinusoidal oscillation of 
which corresponds to a finite excitation energy. When the 
angular momenta of <p(r) and ip(r) are equal, quasiparti- 
cle exchange between them may occur without resulting 
in longitudinal motion of vortices. The finite frequency of 
a mode which emerges due to a spontaneous symmetry 
breaking is typically attributed to the Anderson-Higgs 
mechanism, recently observed in neutral superfluids |48L 
I5U] , Here we observe an opposite effect — the excitation 
frequency of a Tkachenko mode which due to the broken 
rotation symmetry would be expected to be zero but is 
nevertheless observed to be finite under normal circum- 
stances, becomes vanishingly small. 

The Kohn mode or the center-of-mass dipole mode 
with frequency lu± — has been marked in Fig. [3] by 
a triangle. The states with q = 2 — 8 are the lowest ex- 
citation modes in each of the seven Kelvin-Tkachenko 
vortex wave branches present in this system due to the 
seven vortex degrees of freedom [JS]. The two insets 
in Fig. [3] show density isosurface plots of the two zero 
energy states <j>(r) and ip(r) with q — 1 and q = 2, 
respectively. Fig. [2]ja) shows a density isosurface plot 
iV^r)! 2 = |u 2 ( r ) 1 2 = l w 2( r )| 2 viewed along the rota- 
tion axis. Fig. [2jb) shows the phase map M (r) = 
arg(0(r)) — arg(V>(r)). The white markers denote the lo- 
cations of vortex cores in the Majorana boson state ^>(r) 
and the black markers correspond to the vortex phase sin- 
gularities in the condensate ground state. Both 0(r) and 
?A(r) have a vortex phase singularity at the origin, which 
cancel each other in the relative phase map shown in 
Fig. [2jb). Finally, figure [4] shows the probability density 
of the zero energy Majorana boson in the z — plane. It 
can be viewed as a bound state of the effective potential 
created by the combination of the harmonic trap and the 
condensate mean-field. There are six prominent peaks 
which are localized inside the vortex cores of the conden- 



sate. There are also seven vortex phase singularities in 
V>(r). Six of these are located in between the peaks ob- 
served in Fig. [4] and the seventh is located on the z-axis 
in the center of the figure. 

In conclusion, we have found a zero energy quasiparti- 
cle state in the elementary excitation spectrum of a rotat- 
ing Bose-Einstein condensate in a harmonic trap hosting 
seven singly quantized vortices. The physical character 
of this zero mode is that of a Tkachenko vortex wave. 
At zero frequency there is no energy cost associated with 
populating this excitation mode, suggesting the presence 
of a Bose-Einstein condensate of quanta of Tkachenko 
waves. The ground state of this system is doubly degen- 
erate implying a finite entropy of the zero temperature 
ground state. For = wj_ the quantum Hall effect is 
predicted to appear in rotating quantum gases with a 
lowest Landau level and a massively degenerate ground 
state. Here we observe ground state degeneracy already 
at finite uj± — Q. There remains an important question 
regarding the robustness of this zero energy Majorana 
boson mode with respect to quantum and thermal fluctu- 
ations and perturbations in external trapping potential. 

In the Feynman interpretation of antiparticles |51l I52j , 
a Majorana zero mode for which a particle and its an- 
tiparticle are equivalent could be described in terms of 
simultaneous forward and backward time propagation of 
an excitation. Applying this reasoning to the quasiparti- 
cle zero mode in our system and noting that time-reversal 
changes the sense of rotation, we may view the zero en- 
ergy Tkachenko mode as a quasiparticle which comprises 
of simultaneous left and right handed rotation. 

The zero energy Tkachenko mode has the exceptional 
character of corresponding to purely transverse motion 
of the vortices with zero amplitude in the radial com- 
ponent. In contrast, the lowest Tkachenko mode generi- 
cally has a finite excitation frequency consistent with the 
longitudinal mode component which emerges due to its 
coupling to the condensed boson field. This can be com- 
pared with the Anderson-Higgs mechanism in supercon- 
ductors, which is also thought to be responsible for the 
mass generation of elementary particles [151120) . 

I am grateful to Michael Morgan for valuable discus- 
sions. 
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